THE VALUE DISTRIBUTION OF INCOMPLETE GAUSS SUMS 
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Abstract. It is well known that the classical Gauss sum, normalized by the square-root 
number of terms, takes only finitely many values. If one restricts the range of summation to a 
subinterval, a much richer structure emerges. We prove a limit law for the value distribution 
of such incomplete Gauss sums. The limit distribution is given by the distribution of a 
certain family of periodic functions. Our results complement Oskolkov's pointwise bounds 
for incomplete Gauss sums as well as the limit theorems for quadratic Weyl sums (theta sums) 
due to Jurkat and van Home and the second author. 



1. Introduction 

The present paper investigates the asymptotic distribution of the incomplete Gauss sum 

(i-i) g<p(p,q) = XM" ) e i(P h2 )i 

h=0 ^ ^ ' 

where q £ N, p € Z 9 , and e q (x) = e 2nlx / q ; the weight function ip is periodic with period one. 
The case <p = 1 corresponds to the classical Gauss sum. The main example of an incomplete 
Gauss sum in the literature is the case when ip is the characteristic function of a subinterval 
of the unit interval HQ El CGS El QU Q2] . 

It is natural to assume that p and q are coprime, i.e, p£ZJ = {p < q, gcd(p, q) = 1}. Here 
Z* is the multiplicative group of integers mod q. The order of Z* is denoted by 4>(q) (Euler's 
totient function). If p, q are not coprime, say gcd(p, q) = r for some r > 1, we set p' = p/r, 
q' = q/r, and observe that 

(i-2) 9<p(p,q) =sw(pW) 

where 

r-l 

(1.3) (f r (x) = J2if 

k=0 

The case when p, q are not coprime can therefore be reduced to the coprime case. 

Functional equations and pointwise estimates of incomplete Gauss sums have been studied 
extensively by Oskolkov [15] . and the aim of the present paper is to complement his results 
by establishing limit theorems for their value distribution at random argument. 

The existence of a limit distribution of the classical theta sum 

N-l 

(1.4) S N (x) = -= £ <n 2 x), e(x) = e 2 ™, 

* n=0 

for x uniformly distributed in T = M/Z has been proved by Jurkat and van Home [8, 9, 10J (for 
its absolute value) and the second author [13] (for its full distribution in the complex plane) ; 
we refer the reader also to the recent study by Cellarosi [2j. A striking feature of the limit 
distribution of theta sums is that it has a heavy tail: The probability that |iSjv(aO| nas a va l ue 
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greater than R, decays, for large R, as R~ 4 . At rational x, the theta sum of course reduces 
to an incomplete Gauss sum where cp is the characteristic function of an interval, and we will 
see below (Remark [2]) that in this case the limit distribution has compact support — the exact 
opposite of a heavy tail. 

We denote by ^fcez^ e (^ x ) the Fourier series of ip. We will focus for the major part of 
this paper on Gauss sums with differentiable weight functions ip in the space 



(1.5) 



B(T) 



k 2 \pk\ < 00 



and only later extend our results to general Riemann integrable functions, under an additional 
assumption on q. 

The limit distributions of incomplete Gauss sums will be characterized by the following 
random variables: 

• X takes the four values ±1 ± i with equal probability. 

• Y takes the values ±1 with equal probability. 

• Z takes the values 1 ± i with equal probability. 

• G+, Gtp, G~ are random variables given by the Fourier series 

(1.6) G+(x) = ^2p2n e(n 2 x), Gp(x) = ^ (p n e(n 2 x), G~(x) = ^ p n e(n 2 x), 
respectively, with x uniformly distributed on T. 



Remark 1. Note that, for tp G i3(T), the functions in ( |1.6[ ) are differentiable and thus contin- 
uous. If ip satisfies the functional relation p(x) = — x), then its Fourier coefficients are 
related via (p- n = (—l) n ip n , and hence = and = G v . If (p ri + p- n is real-valued for 



all re, then ImG^(-i) = -ImG^(i) and ImG ¥ 



lmG ip (x). Hence the probability 



density describing the distribution of the imaginary part of the random variables G+, G v , G~ 
is symmetric. Furthermore, lmG~(x + \) = ReG^(x), and thus the real and imaginary part 
of G~ have the same distribution. 

We define e a = 1 if a = 1 mod 4, and e a = i if a = 3 mod 4. The symbol A denotes 
convergence in distribution. 

Theorem 1. Fix a subset PcT with boundary of measure zero, and let p 6 B(T). For each 
q £ N, choose p £ Z* n qT> at random with uniform probability. Then, as q —> 00 along an 
appropriate subsequence as specified below, we have: 





q is not a square 


q is a square 


q = mod 4 


1 


r 9i{p,q) 9<e(p,<l)\ d ) , x c+ , 

k Vq 'gi(p,q)J ' v 


( 9i(P,q) g v (p^q)\ d (r7 

I V~q 'mMl^^'^ 


q = 1 mod 2 






e qVq 




q/2 is not a square 


q/2 is a square 


q = 2 mod 4 


fg 1 (2p,q/2) g^(p,q) \ d 


9<p(P,q) d c - 


V e q/2y ft/2 '2 gi (2p, q /2)J 1 


£ q /2V2q v 



The proof of this theorem is given in Section |4j The key ingredients are functional equations 
for incomplete Gauss sums (Section[2]) and estimates on (twisted) Kloosterman sums and Salie 
sums (Section |3]). It is crucial that the exponential sums considered here are quadratic in h. 
The case of higher powers is significantly more difficult, cf. |14j . 
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To illustrate the statement of Theorem [TJ let us consider the distribution of the absolute 
values of incomplete Gauss sums on the positive axis M>o- 

Corollary 1. Under the assumptions of Theorem^ as q —> oo, 



q = mod 4 


1 ^ 


q = 1 mod 2 


Mp,«)I d i 


q = 2 mod 4 


\g v (p,Q)\ d |c _. 



Since for smooth 93 the normalized incomplete Gauss sums q 1 ^ 2 \g v (p, q)\ are bounded 
(cf . Theorem [3] below) , the previous corollary implies convergence of the fcth moment 

(1-7) M k „{q) = ^ |^(p,<z)| fc . 

Corollary 2. Under the assumptions of Theorem^ we have for any k > 



q = mod 4 


lim Mfc - (9) = ) 
9^00 (2g) fc / 2 J 


f\G+(x)\ k dx 

T 


9=1 mod 2 


r M kjip (q) 
hm — 7775 — = / 

q— >oo J 


T 


q = 2 mod 4 


hm Mfc - (9) = ) 
9^00 (2g) fc / 2 y 


f'\G-(x)\ k dx 

T 



The following technical estimate allows us to extend Theorem [TJ to non-smooth <p as long 
as q has a bounded number of divisors d(q). 

Lemma 1. Fix a positive integer N. Then there exists a constant CV > such that, for every 
Riemann integrable function ip : T — )■ C, we have 

(1.8) hmsup < ^ 

d(q)<N 

Lemma [TJ is proved in Section [5j Together with Chebyshev's inequality it implies the 
following extension of Theorem [TJ 

Theorem 2. Fix a subset D C T rai/i boundary of measure zero, and let p : T — >• C be 
Riemann integrable. Then the conclusions of Theorem^ and Corollary^ remain valid for any 
sequence of q — > 00 with a bounded number of divisors. 

The proof of Theorem [2] is supplied in Section [6j 

Remark 2. For ip £ B(T), the functions in ( |1.6[ ) are continuous and bounded. Arkhipov 
and Oskolkov [TJ [15] prove that boundedness (but not continuity) still holds even if ip is the 
characteristic function of a subinterval of T. This implies that the limit distribution has 
compact support and that therefore Corollary [2] remains valid in this case, subject to the 
addtional assumption d(q) < N. 

Remark 3. It seems plausible that the hypothesis on the number of divisors of q can be 
removed from Theorem [2j at least in the case of weights ip of bounded variation — or indeed 
all Riemann integrable functions. 
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-0.2 -0.1 0.0 0.1 0.2 0.3 -0.2 -0.1 0.0 0.1 0.2 0.3 



Figure 1. The histogram on the left shows the value distribution of the real 

part of 9,p< f ,q \ — L where ip is the characteristic function of the interval [0, -4=1, 
1 gi(p,q) y/7 r L ' V7 17 

q = 5012 = 2 2 x 7 x 179 and p uniformly distributed in Z*. There are thus 
4>(q) = 2136 sample points distributed across 40 bins, which means we have on 
average 53.4 values in each bin. The histogram seems consistent with fluctua- 
tions of the order of the square-root of that number, i.e., approximately 14% of 
the height of each bin. The histogram on the right shows the imaginary part 
of the corresponding quantities. The continuous curves represent a numerical 
approximation to the real and imaginary part of the probability density of the 
random variable 



2.5 - 




-0.4 -0.2 0.0 0.2 0.4 0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 



Figure 2. Same as Figure [TJ now for q = 5013 = 3 2 x 557. The number of 
sample points is now <j>(q) = 3336. The continuous curves represent a numerical 
approximation to the real and imaginary part of the probability density of the 
random variable G v {x). 

To illustrate Theorem [2j we have computed numerically the value distribution of the real 
and imaginary parts of incomplete Gauss sums for different values of q, see Figures [l}]3] and 
Section 



2. Functional equations for incomplete Gauss sums 
Legendre's quadratic residue symbol is defined for an odd prime b by 

+1 if b \ a and a is a quadratic residue 



(2.1) 



if b | a 

— 1 if b \ a and a is a quadratic non residue. 



Following Jacobi, we extend the definition to arbitrary odd integers b multiplicatively: Let 
b a positive odd integer with prime factorization J|* =1 p[ l . For a € Z we define (j) = 1 
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Figure 3. The histogram on the left shows the value distribution of the real 
P art °f 2gi(2p'g/2) w * tn tne same V 9 as m Figure jlj and q = 5014 = 2 x 23 x 109, 
p uniformly distributed in Z* , where <j)(q) = 2376. The histogram on the right 
shows the imaginary part of the corresponding quantities. The continuous 
curves represent a numerical approximation to the real and imaginary part of 
the probability density of the random variable G~(x). 



and (f) = Y\i=i(^) ri - The generalized quadratic residue symbol (or Jacobi symbol) (|) is 
characterized by the following properties (cf. [18j): 

(i) (f)=0ifgcd( O) 6)^l. 

(ii) If b is an odd prime, (|) coincides with the ordinary quadratic residue symbol. 

(iii) If b > 0, (j) defines a character modulo b. 

(iv) If a 0, (-) defines a character modulo a divisor of 4a, whose conductor is the 
conductor of Ofv/a) over O. 



(v) (^) 



VI 



sgna. 

1. 

>2 



In particular (|) =1, if gcd(a, b) = 1. 

We assume from now on that gcd(p, q) = 1. 
The classical Gauss sum 

(2-2) 9i(p,q)= e ^h 2 ). 

h mod q 

can be evaluated explicitly in terms of the Jacobi symbol: 

'(l + i)^ 1 ^) VQ if 9^0 mod 4 



(2.3) 



gi{p, q) 



if q = 1 mod 2 
if q = 2 mod 4, 

with e a as defined before Theorem [TJ 

The following theorem is implicit in the papers of Fiedler, Jurkat and Korner [7] and 
Oskolkov [TS] 

Theorem 3. For cp G B(T), 

'<7i(p, <?)<?+(-§) 



(2.4) 



9<p{p,q) 



gi(p,q)G v ( - f) 



8p N 
9/2/ 



if q = mod 4 
if q = 1 mod 2 
if q = 2 mod 4. 



^2 5l (2p,g/2)G- 

(7n t/ie /irsi and second case, x denotes the inverse of x mod g, in £ae i/urd i/ie inverse mod 
a/2j 
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Proof. Since the Fourier series of (p is absolutely convergent, we may assume without loss 
of generality that ip(x) = e(kx) with k G Z fixed. The proof is a then simple exercise in 
completing the square (cf. [7]): 

(i) q = mod 4: For k = 2n even, 

(2.5) E e q (ph 2 + 2nh) = E e g (p(h + np) 2 - p(np) 2 ) = gi(p,q) e q (-pn 2 ). 

ft mod q ft mod q 

For k odd, 

E e q (ph 2 + kh)= e *(p( h + l) +k ( h+ l)) 

^2 h mod <j ft mod <f 

= — E e q (ph 2 + fcfr), 

ft mod <j 

and therefore must be zero. 

(ii) q = l mod 2: For q odd, we may use the inverse of 2 mod q: 

(2.7) ]T e q (ph 2 + kh) = e q (p(h + 2p-k) 2 -p(2l)k) 2 )=g 1 (p,q)e q (-4 1 )k 2 ). 

ft mod <j ft mod q 

(iii) q = 2 mod 4: We deduce claim (iii) from the previous case: Define qo = q/2 and 
po = \{2p - q). Clearly q = 1 mod 2, gcd(p ,<?o) = 1 and | = ^ + |. Then 

(2.8) 

£ + = E 4m 2 + ^ + ^) = E J P0 h 2 + ^^). 

ft mod g ft mod 2qo ft mod 2qo 



Hence, for k odd, A; + go is even and eq. (2.7) yields that the right hand side of d2.8h equals 



(2.9) 



2fll(P0) 90)650 ( - 4p f—^J ) = 2 5i(fo, 9o)e go ( - I6p (k + q ) 2 ) 



^9i(Po,Qo)e qo (- 16p k 



2\ 



2 gi (2p,q/2)e q/2 (-8pk 2 ). 

^2) = (qj2 

Z' 1,2 , (k + q )h\ / - , 2 , (& + <7o)(/i + <7o) 
2^ %)(Po/i + g )= 2^ e 9o(Po(/i + 9o) + ^ 

2^ e go (p /i + = J' 



Furthermore, gi(2p,q/2) = gi(2p,q/2), since (^75) = (^fo)- If & is even, 



ft mod 2(7o v ft mod 2<7n 

(2.10) 



ft mod 2(7o 

This term therefore vanishes. □ 

3. Equidistribution MOD q 

The functional equations of incomplete Gauss sums stated in Theorem [3] lead us to consider 
the joint distribution of | and | on the torus T 2 . The following statement is the second main 
ingredient in the proof of Theorem [T] 

Theorem 4. Let f € C(T 2 ). Then the following convergence holds uniformly in t G Z* as 
g — )• 00: 

(i) For any sequence of q, 

(3-1) :zL£/(?.?W 



9 9/ JT2 



</ 
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(ii) If q = mod 4 is not a square then, for every a E {±1, ±i} ; 

(iii) If q = mod 4 £/ien, /or every a G {±1}, 



p=cr mod 4 

(iv) If q=\ mod 2 is not a square then, for every a G {±1}, 
(3-4) -i- £ f(**)-+l[ f( x )dx. 

(2)= CT 

V q ' 

Remark 4. The statement of Theorem 0] also holds for the test function 
(3.5) f(xi,x 2 ) = Xv(xi)g(,x 2 ), 

where xt> is the characteristic function of a subset PcT with boundary of measure zero and 
g G C(T). This follows from a standard approximation argument. 

Remark 5. The proof of Theorem [4] exploits the classic Weil bounds on (twisted) Kloosterman 
and Salie sums. These directly yield explicit bounds on the rate of convergence in Theorem [4] 
for smooth test functions /. 

The following two lemmas will be helpful in proving Theorem [4j 

Lemma 2. Assume q E N is not a square. Then there exists r G Z such that r = 1 mod 4 
ond(2) = -l. 

Proof. Let (/ = a a 1 1 ...a* s be the prime factorization of q, with ao = 2. Since g is not a square, 
at least one of the U is odd. Suppose now this happens for index j, i.e., tj is odd, for some 
< j < s. 

Fix an integer cq such that, if j = then cq = 5 mod 8 and otherwise Co = 1 mod 8. Then 
{^) = — 1 if i = and otherwise (~) = 1. In both cases of course cq = 1 mod 4. 

Furthermore, let ci, . . . , c s be integers so that for i ^ j, a is a quadratic residue mod aj, 
and for i = j, it is not. Hence (— ) = — 1 if i = j and (— ) = 1 otherwise. 

By the Chinese Remainder Theorem there is an integer r = 1 mod 4 such that (-) = (^) 
and (^r) = (^). Using quadratic reciprocity, we have for i ^ 

since r = 1 mod 4. Using the multiplicativity of the Jacobi symbol we obtain 

□ 



r / V r / [ r J \ r J 



Lemma 3. (i) If q = mod 4 is not a square and a E {±1, ±i}, t/ien 

E 1 = ^(9). 



8 
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(ii) If q = mod 4 and a G {±1}, then 

£ 1 = 



^9 

p=tr mod 4 



(iii) If q = 1 mod 2 is noi a square and a £ {±1}, t/ien 



pez* 

(E)=tr 



Proof, (i) We have 



k Q 



P ez, x feez 4 peZ * KF 

For /c = 0, we have of course, 

(3.9) £l = 0(g). 

pez* 

In the case fc = 1, we need to show that 

(3.10) E*p(1)=°- 

pez* V/V 

By Lemma [2j we can find an r such that (2) = —1 and r = 1 mod 4. Since (') = —1 we have 
r £ and thus for every p £ Z* there is p 6 Z* such that = p. Therefore, 

E = E 

(am ^ mt 

— tr v 



pel 



Since r = 1 mod 4, we have e rp - = ep. Hence, 



(3.12) E% =-E« H 

The case fc = — 1 is the complex conjugate of the case k = 1. For k = 2, 
(3-13) E^) 2 =E^ = 0, 

since = 1 or — 1 if p = 1 or 3 mod 4, respectively. 

(ii) The proof follows from the k = 2 part of the proof for (i), since the assumption that q 
is not a square is not relevant in this case. 

(iii) We have for a = ±1, 



(3-i4) £ 1 = \ £ °- k £ 



9 - F^q 

(H)=cr 
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9 



For k = 0, the sum is obviously equal to 4>{q). The case k = 1 corresponds to the well known 
identity 

(3.15, £(?)=0. 

□ 

Proof of Theorem U We start with the most difficult case (ii). To prove this claim, it 
suffices (in view of Lemma [3] (ii) and Weyl's criterion) to show that 

(3.16) lim-i. Y e ( m P + nt P 

for every fixed (m, n) G 1? \ {(0, 0)}, uniformly in t £ Z* . We have 

— \ / \ k / — 

mp + ntp \ y—^ _ k y—^ k I l\ I m P + 



(3.17) £ ) = E'-*E4.„ 

For A: = the inner sum is the Kloosterman sum 

(3.18) %,nt, g ) = ^e( mP + ^ 

for which we have the classical Weil bound \K(m, nt, q)\ < gcd(m,nt,q) 1 l 2 q l l 2 T(q), see [5]. 
Since m and n are fixed and gcd(t,q) = 1, gcd(m, nt, q) 1 ^ 2 = gcd(m,n,q) 1 / 2 is bounded 
above. Furthermore r(q) <C e q € for any e > 0. Since <f){q) S> £ g 1_e for any e > 0, we see that 
4>{q)~ l K{m,nt,q) tends to zero, uniformly in t 6 Z* as g — > oo, as required. 
The case k = 1 (A; = —1) leads to the twisted Kloosterman sum 

(3.19) s 9 ( m ,„ M ) = £ 



(and the complex conjugate of So(—m, —nt, q)). Here we have the same bound as for Kloost- 
erman sums, \Se(m, nt, q)\ < gcd(m,n,q) 1 ^ 2 q 1 ^ 2 T(q), see [3] (but also the more recent [1]). We 
conclude that the contribution of the k = ±1 term also tends to zero uniformly in t. 
The case k = 2 leads to 

m 2 fmp + ntp\ y-^ (p-l . mp + ntp\ ( q 

(3.20) ^ e p e( j = ^ el — + j = -iK ( m + -, nt, 

and is thus reduced to Kloosterman sums. This proves the case (ii). 
Case (iii) reduces to the same estimates as in case (ii) k = 2. 
Case (iv) is analogous, but here the estimates reduce to bounds on Salie sums 

(3.21) s(m ,*,)_£(E) e (!!£±2»), 



pet 



'•i 



which are the same as the above for the (twisted) Kloosterman sums. 

Case (i) of course follows from the classical Weil bound for Kloosterman sums. □ 



10 EMEK DEMIRCI AKARSU AND JENS MARKLOF 

4. Proof of Theorem CD 

Case la: q = mod 4, q not a square. We need to show that for any bounded continuous 
F : C — > C we have 

In view of Theorem [3j this is equivalent to 

Since and F are continuous, the latter statement follows from Theorem |4] (ii) and subse- 
quent remark, if we choose the test function 

(4.3) f(xi,x 2 ) = xv{xi)F{G+{-x 2 )). 

Case lb: q = mod 4, q is a square. We proceed as in Case lb, and note that the condition 
e p = 1 (e» = i) is equivalent to p = 1 mod 4 (p = —1 mod 4). The statement follows from 
Theorem ffl (hi). 

Case 2a: q = 1 mod 2, g not a square. In this case, the statement to be proved reduces 
(again using Theorem [3]) to 

(") w) E *>OH" f )) - ?//<w- 

which follows from Theorem [4] (iv) with t = 4. 

Case 2b: q = 1 mod 2, q is a square. Analogous to Case 2a, except that we employ 
Theorem [4] (i). 

Case 3a: q = 2 mod 4, g not a square. Following the same strategy as above, we deduce 
that the claim of the theorem is equivalent to 

w> ^ E»©%(-|)^/ T W 



^572) 

As in the proof of Theorem [3] (iii), we substitute g = 2go and p = 2po + qo, i.e., go = <?/2 and 
Po = \(2p — q). Note that this map describes a bijection Z* — > Z* Q . Hence (4.5) is equivalent 
to 



F(G-(x))dx, 



u 10 



which is again implied by Theorem [4] (iv) . 

Case 3b: q = 2 mod 4, q is a square. Analogous to Case 3a, except that we use Theorem 
|(i). ' ' □ 



1. 



5. Mean-square estimates 
The key step in the proof of Theorem [2] is the estimate on the mean-square given in Lemma 
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Proof of Lemma [7} We have 
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(5.1) 



< 



E 



h\ fti 
(pi - )(p 



h 2 =h' 2 mod q 



Take ip £ S(T) such that |^(a;)| < ^0*0 for all x £ T. Then 



h,h'eZ q 
h 2 =h' 2 mod 5 



(5.2) 



r^Ew^)i 2 

J^)E E l^(P,9/r)| S 



where (recall (1.2)) 
(5.3) 



^(*) = X>(^). 

fc=0 v 7 



q/r 



The Fourier series of this function is 

(5.4) ip r (x) = r} j ip rn ei 



nx 



net 



where ip^ are the Fourier coefficients of ip. We have thus shown that 



(5 ' 5 > **M*ji>jk)£*(r) 



.1/ 



(v=t) ( q 



< q V 1 M (P=T) f ? 
- 12310(g) ^ r 2 ^ 



V r 



By Corollary [2j for every fixed r, 

(5.6) lrm -M 2 ^ r (-^j < 2r ($ \ 2 + $ rn + ^ rn \A < 2r (\^\ 2 + \$ n + $- n 

The right hand side is bounded by 4r||^||2. The convergence is uniform in r, if we assume ip 
has a finite Fourier series. In this case, we therefore have 



(5.7) 



lim sup 



M 2tV (q) < C N 



q— >oo q 
d(q)<N 



Since ip is Riemann integrable, given any e > 0, there exist ip with finite Fourier series such 
that (a) \<p\ < ip (as required in (5.2)) and (b) \\\p\ —ip\\2 < £■ This proves that the right hand 
side of (5.7) is arbitrarily close to |MI|- ^ 



6. Proof of Theorem [2] 

The following lemma says that the sequence probability measures defined by the value 
distribution of incomplete Gauss sums is tight. By the Helly-Prokhorov theorem, this means 
that the sequence is relatively compact, i.e., every sequence contains a convergent subsequence. 



12 EMEK DEMIRCI AKARSU AND JENS MARKLOF 

Lemma 4. Fix N . For every e > there exists K e > such that 



(6.1) limsup-^freZg* -.q-^lg v (p,q)\ >K e }\ < e\\<p§ 

q^oo (p(q) 
d(q)<N 

for any Riemann integrable ip : T — > C. 
Proof. We have, by Chebyshev's inequality 

(6-2) ^ ) \{P^q^- 1 %AP^)\>K}\<^. 

The claim now follows from Lemma [TJ □ 



Chebyshev's inequality (6.2) also implies the following 



Lemma 5. Fix N. Let <p : T — >■ C be Riemann integrable. Then, for every e > 0, 5 > there 
exists if) G B(T) and such that 

(6.3) limsup — r^l-fa e Z* : ^ 1/2 |£^(p, <?) - &/>(p> ?)| > #}| < e. 

q^oo (j)(q) 
d(q)<N 



Proof. This follows immediately from (6.2); note that g<p(p,q) — gip(p,q) = g v -ip(p,q) and 



9? — ^ is Riemann integrable. □ 

We now turn to the proof of Theorem[2} We restrict ourselves to Case la where q = mod 4; 
the other cases are analogous. The relative compactness implied by Lemma [4] can be stated 
as follows. Any sequence of q — > oo with d(q) < N contains a subsequence {qj} with the 
property: there is a probability measure v on {±1 ± i} x C such that for any a E {±1 ± i} 
and any bounded continuous function F : C — > C we have 

(6.4) lim — \-- Y f( 9 ^ q A = f F(z)u(a,dz). 

The probability measure f may depend on the choice of subsequence, <p and on P. 

Let us now show that for every F E C§°(C) (infinitely differentiable and of compact support) 
the limit 

(6-5) I V (F):= lim V F^*^ 



exists. Iip(F) must then be equal to the right hand side of (6.4), which in fact means that u 
is unique and the full sequence of q converges. 

To prove the existence of I V {F), note first of all that since F E C§°(C) we have \F(w) — 
F( z )\ < Cmin{l, \w — z\} for some constant C > 0. Therefore, for tp, 5, e as in Lemma[5j we 
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have 



W(9) 



E 



/ g^(g) g) 

\5i(P>g) 



(6.6) 



< 



< 



C 



E min ] 1 > 



|2>|0(<Z) 

c 



^ min i 1, 



3l(p,g) 01 (p>?) 



9<pb,g) 9^ip,Q) 



0l(p>g) 5i(P,g) 



< M (2 1/2 ^ + 6). 



Since the limit I^(F) exists by Theorem [TJ the sequence 



(6.7) 



l^l<Mg) 



E 



gy-(p, g) 
9i(p,q) 



defines a Cauchy sequence. Using this fact, the bound (6.6) and the triangle inequality, we 
see that 



1 



\v\<Kq) 



E * 



( 9y{p, g) 
V5l(P) g) 



is a Cauchy sequence, too, and hence Iip(F) exists. As mentioned earlier, this means that 
Itp(F) must then be equal to the right hand side of ( |6.4| ), which in fact means that v is unique 
and the full sequence of q converges for every bounded continuous F. 

The bound (6.6) furthermore implies that I^(F) — > I<p(F), as ip — > ip in L 2 (T). This 
completes the proof of Theorem [2j 



7. Numerics 

The computations used in Figures [l]-[3] were carried out with Mathematica. We encoded the 
real and imaginary part of the incomplete Gauss sum — ^ (where cp is the characteristic 

function of the interval [0, ^] C [0, 1]) as 

ReGauss[p_, q_, T_] := 
If[GCD[p, q] == 1, 
Re[Sum[Exp[2*Pi*I*rT2*p/q] , {h, 1, T>] / 
Sum[Exp[2*Pi*I*h~2*p/q] , {h, 1, q>]] - T/q, Infinity] 
ImGauss[p_, q_, T_] := 
If[GCD[p, q] == 1, 
Im[Sum[Exp[2*Pi*I*rT2*p/q] , {h, 1, T>] / 
Sum[Exp[2*Pi*I*h~2*p/q] , {h, 1, q}]], Infinity] 

and formed a table comprising the values for all integers p < q. Whenever gcd(p, q) ^ 1 the 
value oo is assigned, which is ignored by Mathematica's Histogram command. 

The probability density of real/imaginary part of G^ and G v in Figures 111 and El was plotted 
via the SmoothHistogram command, where we truncated the Fourier series G^(x) and G ip (x) 
at n = 4000 and sampled x at 300,000 random points in [0, 1]. As the distribution of real 
and imaginary part of G~ are the same, we only computed ImG^J in Figure [3j truncated at 
n = 5000 and with 500,000 sample points. 
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